function U1 = U1(c,gamma,cf)


 if gamma==1
U1=cf*log(c)+0*gamma;
 else
     U1=cf*c^(1-gamma)/(1-gamma);
 end     


end